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NEW CONVERGENCE CRITERIA FOR THE VORTEX-LATTICE MODELS 
OF THE LEADING-EDGE SEPARATION* 

Osama A. Kandil, Dean T. Mook, and All H. Nayfeh 
Virginia Polytechnic Institute and State U Iversity 

SUMMARY 

The convergence criterion for the vortex-lattice technique which deals with 
delta wings exhibiting significant leading-edge separation has two re luirements. 
First, the wake must converge to a force-free position for any given number of 
discrete vortex elements. Second, the distributed loads must converge as the 
number of elements increases. Replacing the vortex sheets representing the 
wakes by a system of discrete vortex lines whose positions are determined as 
part of the solution (first requirement), one finds that the total loads com- 
puted agree very well with experimental data. But the predicted pressure dis- 
tributions have some irregularities which are the result of discrete vortex 
lines coming close to the lifting surface. Here it is shown that one can 
eliminate these irregularities and predict pressure distributions which agree 
fairly well with experimental data (which show .some irregularities of their 
own) by replacing the system of discrete vortex lines with a single concentra- 
ted core. This core has a circulation equal to the algebraic sur. of the cir- 
culations around the discrete lines and is located at the centroid of these 
lines. Moreover, the second requirement is replaced by the requirement that 
the position and strength of the core converge as the number of elements in- 
creases. Because the calculation of the position and strength of the core is 
much less involved than the calculation of the loads, this approach has the ad- 
ditional desirable feature of requiring less computational time. 

INTRODUCTION 

A characteristic feature of the flow over wings having highly swept, sharp 
leading edges is the formation of vortices above suction uides in the vicinity 
of the leading edges. These vortices roll up in a conical-like spiral with a 
concentrated core. This vortex spiral grows in size and strength as it approaches 
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Che trailing edge. Below the angle of stall (i.e., the angle at which vortex 
bursting occurs), the effect of the leading-edge separation Is to increase the 
velocity on the suction side of the lifting surface by adding a strong cross- 
flow component and hence to increase the aerodynamic loads. Experimental in- 
vestigations, such as those described in refs. 1-7, and numerical investigations, 
such as those described in refs. 8 and 9, confirm these conclusions. Detailed 
descriptions of the flow field are als' given in these references. 

Though the experimental results show that the pressure distribution is some- 
what influenced by the character (laminar or turbulent) of the boundary layer, 
the lift and pitching-moment coefficients are independent of the Reynolds number. 
Thus, one expects an invlscld model of the flow to predict the total loads more 
reliably than tne distributed loads. And such has turned out to be the case. 

The early attempts to develop an invlscld model of these flows were based 
on assumptions of conical flow and/or slender-body theory. These assumptions 
are apparently discredited by the experimental observations; more discussion is 
given in references 8 and 9. 

Subsequent attempts to develop inviscid models did not use these assumptions. 
Rehbach (ref. 10) developed a vortex-lattice technique which progressively short- 
ens the leading edge of a rectangular wing until a delta wing is formed. Associ- 
ated with this method are questions concerning the second convergence requirement 
and the undesirable feature of long computation times; mere discussion is given 
in refs. 8 and 9. Weber et al (ref. 11) developed a technique which uses piece- 
wise continuous, quadratic, doublet-sheet distributions. The second requirement 
apparently was not considered, and no results showing the computed shapes of the 
rol.led-up wakes were given. 

In a related effort, Suciu and Morino (ref. 12) developed a technique for 
modelling the region adjoining the trailing edge. However, numerical experiments, 
as described in refs. 8 aid 9, show that the wake adjoining the trailing edge has 
very little influence on the aerodynamic loads. 

In the present paper, the question of the convergence of the centro idal 
line as the number of discrete elements increases is considered; and this line, 
instead of the system of discrete vortex lines, is used to compute the pressure. 
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aspect ratio 
wing semi-span 

pitching-moment coefficient about the z-axis 
normal-force coefficient 
local cross-normal-force coefficient 
pressure coefficient 
root chord 

local semi-span/wing semi-span 
thickness ratio 
free-stream velocity 

body-fixed axes (z-axis is in spanwise direction) 

dimensionless chord station 

dimensionless spanwise station 

angle of attack 

dimensionless circulation 

number of elemental areas of the lattice 


FORMULATION OF THE PROBLEM 

The perturbation velocity ootential, <J>, of the inviscid* irrotational, in- 
compressible flow past a wing is governed by Laplace’s equation and satisfies 
the following boundary conditions: (a) the no-penetration condition on the wing 
surface given by 

(U 4* V<t>) • n * 0 on S(r) = 0 

00 

where n is the unit normal to the wing surface S, (b) the no-pressure discontin- 
uity condition across the wakes emanating from the leading and trailing edges of 
the wing given by 

Ap = 0 across w(r) = 0 

where Ap is the pressure jump across the wake surface w. This surface is an un- 
known of the problem and must be obtained as part of the solution, (c) the Kutta 
condition which requires that no-pressure jump exists across the wing surface 
along the leading and trailing edges where the wake surface is emanating, and 
(d) the disturbance velocity, V<j>, diminishes far from the wing surface, S, and 
the wake surface, w. We note that the problem is nonlinear due to boundary con- 
dition (b). 


287 



r 


i— i t 1 


DESCRIPTION OF THE METHOD OF SOLUTION 


The solution of the problem posed above is constructed by modelling the 
lifting surface with a bound-vortex lattice and the wake with a system of dis- 
crete, nonintersecting vortex lines. Each vortex line in the wake is composed 
of a series of short, straight segments and one final semi-infinite segment. The 
unknowns here are the circulations around the vortex segments and the positions 
of the finite segments in the wake. 

The disturbance velocity field produced by this model of the wing and wake 
is calculated according to the Biot-Savart law; thus, everywhere, except on the 
wing and the wake, Laplace's equation is satisfied. Moreover, the disturbance 
created by the wing dies out far from the wing and wake, boundary condition (d) . 

Associated with each elemental area of the lattice and with each finite seg- 
ment in the wake is a control point. The lattice is arranged so that vortex seg- 
ments leaving the sharp edges do so at right angles to the edge. Moreover, con- 
trol points are placed between the last edgewise vortex segment of each row and 
column and the edge itself. This arrangement partially satisfies the Kutta con- 
dition (c) . 

The circulations and the positions of the finite segments in the wake are 
obtained by simultaneously requiring the normal component of velocity to vanish 
at the control points of the elemental areas of the lattice and the finite seg- 
ments in the wake to be parallel to the velocities at their own control points. 
Boundary conditions (a) and (b) and the Kutta condition (c) are then satisfied, 
and the problem is solved for this lattice (i.e., the first requirement of the 
convergence criterion is met). More details, especially those regarding the 
iterative, procedure used to effect the last step, are given in refs. 8 and 9. 

Instead of calculating the aerodynamic loads and testing the second re- 
quirement of convergence at this point, with the present procedure we calculate 
the centroidal line of the system of free-vortex lines representing the leading- 
edge wake. To construct this centroidal line, we consider a series of cross- 
flow planes. Proceeding from the vertex toward the trailing edge, we calculate 
the centroids of the vortex lines penetrating these planes according to 
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where r , is taken to be the circulation around the centroidal line between the 
cl 

j*"* 1 and (j+l)^ planes and it given by 


v j 

r = F r. 

cJ i=i 1 

is the position of the centroid in the j 11 * 1 plane, r_ is the position of the 

intersection of the i ^ vortex line with the plane, F. is the circulation 

t h ^ t h 

around the i line, and n_. is the number of lines penetrating the j plane. 

More lines penetrate the planes near the trailing edge than those near the ver- 
tex; thus, T c j increases toward the trailing edge. 

Now, the number of elements is increased and new centroidal lines are cal- 
culated until the changes in r ^ and r^ fall within prescribed tolerances. At 
this point, the second requirement of the convergence criterion is met. 

Only when both requirements are met, do we calculate the aerodynamic loads 
and pressure distribution. The details for calculating the loads are given in 
ref. 8. The numerical results below show, for the examples being considered at 
least, that the centroidal lines converge to a position which is very close to 
the experimentally determined position of the core, that the total loads agree 
very well with the experimental data, and that the pressure distributions agree 
fairly well with the experimental data. 


NUMERICAL EXAMPLES 


Figures 1 and 2 show the actual calculated positions of the free-vortex 
lines for two delta wings. The plan view also shows the bound-vortex lattice. 

And the three dimensional view shows the free-vortex lines, their centroidal 
line, and the trace of the spiral vortex sheet. All the following results are 
associated with these two wings. 

Figures 3 and 4 show the convergence of the circulations around the cen- 
troidal line. And figures 5-8 show the convergence of the position of the cen- 
troidal line and the close agreement between the position of the calculated line 
and the experimentally determined position of the vortex core. 

Figures 9 and 10 show the convergence of the total loads calculated by using 
the system of discrete lines and by using the centroidal line as a function of 
the number of elements (nel) . These results are compared with those obtained by 
the leading-edge-suction analogy (ref. 13) and with experimental data. We note 
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that there is a considerable difference in the experimentally determined normal- 
force coefficients in figure 10. 

Figures 11 and 12 show comparisons of the predicted cross-load coefficients 
and experimental data. And figures 13-16 show comparisons of the predicted 
pressure distributions at several chordwise stations with those obtained by 
another method (ref. 11) and with experimental data. We note that the shape and 
size of the suction peak on the upper surface under the vortex differ from one 
experiment to another, depending on how thick the wing is and on whether the 
boundary layer is laminar or turbulent (refs. 5, 6, and 14). 

CONCLUDING REMARKS 

The second requirement of convergence is based on the centroidal line of 
the free-vortex lines representing the wake. Using this requirement greatly 
reduces the computational time. The position of the centroidal line compares 
very well with that of the vortex core. The centroidal line can also be used to 
calculate the total and distributed aerodynamic loads with good accuracy. This 
results in more reduction in the computational time and smoothing of the peaks 
produced by using many discrete lines. 
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Figure 3.- Variation of circulation along the vortex core 
of a delta wing with fo - 1. a = 15°. 
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Figure 4.- Variation of circulation along the vortex core 
of a delta wing with & * 1.46. a * 14° . 
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Figure 5.- Spanwise position of vortex-core path 
on a delta wing with 5L * 1 . ot ® 15°. 
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Figure 6.- Height of vertex-core path on a 
delta wing with &. * 1. a ■ 15°. 
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Figure 9.- Normal- force and pitching-moment 
coefficients of a delta wing with Si « 1. 
a - 15°. 
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Figure 10.- Normal-force and pitching-moment 
coefficients of a delta wing with A * 1.46. 
a - 14°. 
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Figure 11.- Longitudinal distribution of cross load 
for a delta wing with JR * 1. 12x12 lattice; 

a * 15°. 
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Figure 12.- Longitudinal distribution of cross load 
for a delta wing with JR - 1.46. 12x12 lattice; 

a - 1 
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Figure 13.- Surface pressure distribution of a delta 

wing at — - 0.7. /R - 1; a - 15°. 
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Figure 14.- Surface pressure distribution of a delta 
wing at — ■ 0.82. JR * 1: a » 15°. 
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Figure 15.- Surface pressure distribution of a delta 

wing at — = 0.67. At = 1.46; a * 14°. 
c r 
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Figure 16.- Surface pressure distribution of a delta 

wing at — ■ 0.83. At n 1.46; a ■ 14°. 
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